Probing Matter Radii of Neutron-Rich Nuclei by Antiproton Scattering 
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We propose to use antiprotons to investigate the sizes of stable and neutron-rich exotic nuclei 
by measurements of the pA absorption cross section along isotopic chains in inverse kinematics. 
The expected effects are studied theoretically in a microscopic model. The pU optical potentials are 
obtained by folding free space pN scattering amplitudes with HFB ground state densities and solving 
the scattering equations by direct integration. The mass dependence of absorption cross sections 
is found to follow closely the nuclear root-mean-square radii. The total absorption cross section is 
shown to be a superposition of cross sections describing partial absorption on neutrons and protons, 
respectively. Thus measuring the differential cross sections for absorption on neutrons and protons 
will give information on their respective distributions. In neutron-rich nuclei the outer neutron 
■ layer shields the absorption on the protons giving access to investigations of antiproton-neutron 

D ■ interactions in matter. 
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I. INTRODUCTION 

in 

\Q . To understand the evolution of nuclear sizes and shapes from the bottom to the edges of the valley of /3— stability 
is one of the central questions of modern nuclear structure physics. Experimentally, a variety of efforts is undertaken 
or in planning to investigate properties of neutron-rich exotic nuclei ranging from high-energy breakup reactions 1] 
and elastic proton scattering to charge-exchange reactions Q and low-energy transfer reactions 00. Here, we 
propose to use antiprotons to probe nuclear sizes of stable and unstable nuclei in a systematic way. 
— . ■ Investigations of radii and density distributions of stable nuclei by antiprotons are by itself a long discussed and 
\ applied method, e.g. at BNL 0, [?f or LEAR [1,0, El El- Hitherto, the experiments have been performed with 
I . secondary antiproton beams on a variety of stable nuclei in fixed target geometry and standard kinematics, actually 
using stopped antiprotons. Obviously, another approach must be applied if antiprotons should be used for reactions 
on short-lived isotopes. These nuclei by themselves are available also only as secondary beams, produced either 
by fragmentation or isotope separation on line. A solution dissolving these conflic ting conditions is to perform the 
\ measurements in colliding beam geometry. Such a setup was recently proposed in [l^ . It will become feasible with 
the meanwhile approved FAIR facility a GSI ^4|. The FAIR plans include as central components p production and 
accumulation facilities. For nuclear structure research the New Experimental Storage Ring (NESR) for short-lived 
nuclides and an intersecting electron accelerator for colliding beam experiments will be available El • As discussed in 
the Antiproton-Ion-Collider (AIC) proposal 0], with moderate modifications on an acceptable level the NESR/e - 
setup can be converted into a pA collider facility. Details of the experimental setup and procedures for the measurement 
of antiproton absorption cross sections using Schottky noise frequency spectroscopy of the coasting reaction products 
in the NESR are found in the NUSTAR letter of intent and proposal for FAIR [ia[lr|. 

The physics of pA interactions is by itself an interesting topic for nuclear and hadron physics. It was brought to 
the attention of nuclear physics originally by Duerr and Teller El . The usefulness of p + A scattering for studies of 
the neutron skin in stable medium and heavy nuclei was realized some years later in an experiment at BNL In 
the aftermath, a variety of experiments have been performed, e.g. [cl I9L ITol 111] ]. As for the elementary antiproton 
nucleon (pN) vertex the pA reactions are dominated by processes in which the incoming antiproton annihilates on a 
target nucleon into a variety of particles producing typically and preferentially a shower of pions as the final result. 
Here, we utilize the strong pA absorption for nuclear structure investigations without paying particular attention to 
the hadron physics aspects, some of which, however, can - and in another context will - be studied with the facility 
under discussion, e.g. with the PANDA detector, also being part of the FAIR proposal. 

For our purpose we are satisfied to know that the elementary annihilation processes result finally in a strong 
suppression of the incoming flux. In the elementary elastic scattering amplitude fpN the annihilation channels appear 
as a strong imaginary part while only a small, almost vanishing, real part is observed. Correspondingly, in potential 
models for inclusive pA elastic scattering the annihilation channels are typically accounted for globally by a strongly 
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absorptive optical potential, as e.g. in [20|,|2lJ,|22j,|23j. In fact, the total cross section is about twice the elastic cross 
section |24| which underlines the importance of the annihilation channels. 

The new aspect of the present paper is to explore theoretically the perspectives of antiproton scattering on stable 
and rare nuclei and the use of such reactions for nuclear structure investigations. In particular, we consider elastic 
pA scattering and study to what degree the total reaction or absorption cross section provides information on the 
size and shape of the target nucleus. The theoretical approach is described in sect. [HI As a typical and interesting 
case absorption cross sections for the Ni isotopic chain are discussed in sect III Bl The paper closes with a summary 
in sect lllll 



II. MICROSCOPIC DESCRIPTION OF ANTIPROTON ABSORPTION ON NUCLEI 

A. Optical Model Approach for pA Elastic Scattering 

The basic ingredients for a microscopic description of pA scattering are obviously the elementary antiproton- nucleon 
interactions. Including a spin-independent amplitude A{-\fs, t) and a spin-orbit part C(y/s,t) the free space antiproton 
scattering amplitude is taken to be 

fpN(V^,q) = MV^,t) + c(^,t)^-a ¥ -((ixk) , (l) 

where k is the momentum of the incoming antiproton, s the Mandelstam center of mass energy and q 2 = —t is 
the momentum transfer. For kinetic energies Ti a b ^ 1500 MeV the data clearly show that the pN interactions are 
dominated by annihilation processes. By means of the optical theorem an effective T-matrix is defined 

tp N (T Lab , q 2 ) = -Jf^ a pN(T L ab) (1 - ie) F pN {q 2 ) (2) 

where Ttab = ELab{k) — M is the kinetic energy of the incoming antiproton with momentum k and the total pN 
cross section (JpN{TLab)- The ratio of the real and the imaginary part is denoted by e which lies between \ ' ■ ■ \ in 
the energy range considered here. As in [^J we use a t-channel form factor of gaussian shape, F(q 2 ) = er 13 « . A 
closer inspection shows that by eq|21ipAr is assumed to be separable into an on-shcll strength factor and an off-shell 
form factor. 

The coordinate space pA optical potential is then obtained in the impulse approximation and after a Fourier 
transformation, 

U opt (r)= J2 J ^PN(q)tpN(TLab,q 2 )e 1 ^ , (3) 

N—p.n 

a strongly absorptive complex optical potential U opt — V + iW with real and imaginary parts V and W, respectively, 

is obtained. Obviously, we have U opt = llj^l + U^}. 

In practice, we use the momentum space amplitudes, including also a complex spin-orbit part, derived in [2{| 
according to the Paris model . As in the elementary pN interaction the potentials are dominated by strong 
imaginary parts. We also include the (long range) Coulomb potential U c which is calculated microscopically by 
folding the elementary pN Coulomb amplitude with the HFB charge density p c . The elementary proton and neutron 
electric form factors G^' p \q 2 ) are taken into account. The pA spin-orbit optical potentials U so are calculated with 
the HFB spin densities. Although the resulting U so are rather weak and, in particular, do not contribute significantly 
to the absorption they are taken into account in the calculations for completeness. 

We use a non-relativistic approach and describe pA scattering by the optical model Schroedinger equation 

h 2 — * \ 
-^rV 2 + U pt-T cm \ *L +) (k,r) =0 . (4) 

Here, ^i + ^(k, r) is the optical model wave function for the incoming proton in spin state m a = ±^ and with 
asymptotically outgoing spherical waves. The antiproton-nucleus center-of-mass system with reduced mass m is used. 

From the solutions of eq^Jwe determine the partial wave S-matrix elements Sgj |28| where for an absorptive system 
\Stj\ < 1. The continuity equation for the optical model wave function <I>( + ) provides us with the exact relation 

*<*• = ^E|ti( 1 -|%| 2 ) = -||2^TT^(/ ^ (+)t (k,r)S([/ opt )(r)^ + )(k,r)) , (5) 
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with averaging over the p spin projections. The second equality provides us with a particular useful relation: From 
eqOHwe know ^s(U op t) = W p + W n which immediately implies additivity for a a b s = a^ s + cr^ s and the partial cross 
section for annihilation on q — p, n, respectively, are given by 



( g ) _ 2m 1 
°° 6s ~ ~ kh 2 2s + 1 



tr s (^J d 3 r\^(k,r)\ 2 W q (r^j . (6) 
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FIG. 1: Root-mean-square (rms) radii for 48 ~ 88 Ni-isotopes. HFB results for protons (squares), 
neutrons (circles) and the total (isoscalar) density (stars) are shown. The mass range corresponds 
to the isotopes found in the HFB calculations to be particle-stable. The lines are to guide the eye. 
Note the change from proton skins at A < 58 to neutron skins at larger mass numbers. 



For our aim eqEJis of particular interest because the leading order proton and neutron contributions to a a bs are being 
displayed explicitly. Hence, we have succeeded to represent the annihilation on the target protons and neutrons in a 
formally separable form allowing to study theoretically the corresponding partial absorption cross sections separately. 
However, we note that in eqElthe distorted waves are determined by U opt in a non-perturbative manner and, as such, 
include the whole multitude of effects from W P;n . 

The ground state densities entering into the optical potentials are taken from non-relativistic HFB calculations 
using a G-Matrix interaction as described in |29| . The known deficiencies of non-relativistic Brueckner-calculations 
are overcome by introducing additional density dependent contributions pflllSlj and adjusting the parameters to the 
variational results for infinite nuclear matter of the Urbana group [3(| • An effective density dependent zero-range 
pairing interaction is derived from the SE (S=0,T=1) in-medium NN-interaction |3l| . The system of coupled state 
dependent HFB gap equations is solved self-consistently. HFB results for rms-radii of the proton and neutron ground 
state densities in the Ni isotopes are displayed in Fig^ We also find a very satisfactory overall agreement of the 
theoretical and the measured binding energies on the level of 5%, see e.g. [3ll 132^ 



B. Absorption Cross Sections for p Annihilation on the Ni-Isotopes 



While for most pA investigations the strong absorption is an unwanted and disturbing effect, we take advantage 
of this property: Because of the strong absorption the product of scattering waves in eq|S] gains strength mainly 
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in the nuclear surface thus favoring investigations of nuclear skin configurations. As seen from eq^this behavior, 
common for both the proton and neutron partial cross sections, is modulated by the specific shape of W p , n leading 
to a pronounced weighting of the nuclear surface region. Thus, measurements of the pA absorption cross section will 
provide us with data proving the existence of nuclear skins and recording their properties like their composition and 
mass dependence. 

The total absorption cross sections for Ni-Isotopes in the mass range 48<A<88 are shown in FigEl Results for 
antiproton incident energies from T^f, =50 MeV to TLab =400 MeV are displayed. With increasing energy the 
magnitude of a a \, s decreases at a rate slowing down when approaching the highest shown energy, Ti ao — 400 MeV. 
A remarkable feature common to all energies is the steady increase with mass number, reflecting the growth of the 
neutron skin with the neutron excess. 
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FIG. 2: Absorption cross sections for antiproton annihilation on the 48_88 Ni-isotopes at various 
antiproton incident energies. The lines are to guide the eye. 



In Fig. 3 the absorption cross sections at Ti a h=200 MeV and T£ a b=400 MeV are compared to the rms-radii of the 
nuclear ground state number density distributions. At both energies the nuclear rms-radii were normalized arbitrarily 
to the cross sections at 58 Ni expressed by the scaling law cr afcs ~ a (r 2 ). From Fig. 3 we note an impressive agreement 
for (Tabs and the rms radii, both increasing with almost the same slope as functions of the target mass. This behavior 
is found at all energies considered here, thus indicating an almost universal scaling law . The normalization constants 
ceo are slowly varying with the incident energy, decreasing by about 5% from T^j, = 200 MeV to T^ab = 400 MeV. 
Their origin and physical content will be discussed below. We conclude that antiproton annihilation is a well suited 
probe for investigating nuclear sizes and shapes. 

It is instructive to inspect pA scattering along a trajectory through the target nucleus by means of the eikonal 
approximation [28j . The partial absorption cross sections for q = p, n at impact parameter b are given by 

rj — /> + 00 

with the eikonal absorption kernel 

= ^^tr s {Tl^(z,b)} = -e- 2x '' (z ' b) W q (z,b) , (8) 
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FIG. 3: Absorption cross sections (filled squares) for antiproton annihilation on Ni-isotopes at Tiab = 200 MeV and Tiab = 
400 MeV, respectively, are compared to the rms-radii of the nuclear matter densities (filled circles). In both cases, the rms- radii 
are arbitrarily normalized by a factor Qo to a a bs for 58 Ni. Also shown are the partial cross sections for absorption on the target 
neutrons (stars) and protons (diamonds). 



including a trace over spin projections. The eikonal integral 

Xq ^ h) = ~k^ J_ d ^ W ^ b) (9) 

describes the depletion of the incoming antiproton flux along the trajectory through the target nucleus. These relations 
show that n(z, b) is a measure for the degree of absorption at a given point in the (z, b) scattering plane. Overall, 
for impact parameters b less than the target radius Ra the incoming flux is strongly absorbed, giving the target the 
characteristics of a black disk. 

In Fig. 4, results for 11(2, b) are displayed at Tiab = 300 MeV for a variety of Ni nuclei and an impact parameter 
b = 4 fm which roughly corresponds to grazing collisions over the mass region considered in Fig. 4. We find an 
interesting mass dependence: In the neutron-deficient 48 Ni the absorption is takes place dominantly in the proton 
skin and with increasing neutron number the process changes gradually to absorption on the neutrons with a break- 
even point around 58 Ni. In the extremely neutron-rich isotopes beyond 78 Ni the antiprotons are absorbed preferentially 
in the neutron skin, i.e. in the outer layers of the nuclear density distribution. As a consequence, in the neutron- 
rich isotopes the more tightly bound and spatially more confined protons are effectively screened from the incoming 
p flux by the neutron skin. This implies that the antiprotons are preferentially annihilated in a region of almost 
pure neutron matter. The screening effect also explains the decrease of the partial proton and the increase of the 
partial neutron optical model cross sections, respectively, seen in Fig. 3. In fact, the eikonal approach leads to a very 
satisfying agreement with the exact optical model results even on the quantitative level. Above Jz, a (, ~ 100 MeV both 
descriptions agree by better than 10% and, as to be expected, the deviations diminish with increasing incident energy. 
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FIG. 4: Absorption of the incoming antiproton flux along a trajectory with impact parameter b = 4 fm. The eikonal absorption 
kernel, eq|5] is displayed for annihilation of antiprotons with an incident energy of Tt a t = 300 MeV on 48 . 58 > 68 . 78 Ni The p beam 
is incoming from the left, i.e. from z = — oo. Note the reversed behavior of the p absorption on protons and neutrons with 
increasing neutron excess. 



Assuming a gaussian form factor for W(r) the eikonal cross section, eq[7| can be calculated in closed form, resulting 
in a a bs(A,T cm ) = f(A, T cm )-n: (r 2 ) a- The eikonal scaling function 

f(A, T cm ) = 7 + log T cm )) + Ei(£(A, T cm )) (10) 

describes the deviation from the black disk limit a a bs — * ^{r 2 )A- 7 = 0.5772 • • • denotes Euler's constant, Ei(x) is 
the exponential integral [34| and 

&A,T m ) = ^'M% NiT ° m))po . (11) 

cm 

is given in terms of physical quantities, including the pN T-matrix, the central nuclear density po and £ g = kRA 
corresponding to the pA grazing angular momentum. For the physically relevant values of £ the Ei{^) term can 
safely be neglected. These relations enable us to interpret in physical terms the scaling factor introduced above. By 
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comparison we find ao ~ a,Q — ^f(A — 58, T cm ). Indeed, approximating the microscopic W(r) by a Gaussian of the 
same volume integral and rms-radius the analytic eikonal results agree with the full optical model ao convincingly 
well, e.g. ao/ao — 0.92 ■ • • 0.94 for TL a b = 100 ■ • • 400 MeV and unity is approached a higher energies. The functional 
structure of the eikonal scaling factor, ea llOl also explains the apparent convergence of the absorption cross sections 
to limiting asymptotic values as indicated in Fig|2]for the larger energies. Thus, f{A 1 T cm ) and correspondingly ao 
contain valuable information on the dynamics of the pA system. 



III. CONCLUSIONS AND OUTLOOK 

Antiproton-nucleus scattering was investigated theoretically as a tool for measuring neutron skins in exotic nuclei. 
The pA absorption cross sections were calculated with microscopic optical potentials derived from the elementary 
free space pN scattering amplitudes and HFB density distributions. The results point to promising perspectives in 
at least two directions, namely 

• showing that antiprotons are a perfect probe for investigations of neutron (and proton) skins; 

• indicating a new method for investigating antiproton-neutron interactions. 

The latter point is of particular interest because the data base on pn interaction is very scarce. It also adds a new 
aspect to the experiment proposed in 0, : Scanning over an isotopic chain will allow to record the gradual change 
from pp to pn annihilation, hence providing a new access to the p&ar-nucleon interactions. 

Our results show that simple absorption cross section measurements will already provide valuable new information 
on the densities of the target nuclei, most likely of a higher accuracy than achievable by other methods. More involved 
measurements seem to be feasible from which more detailed information on the configurations of the target nuclei can 
be obtained. 

A next step will be to consider more exclusive experiments. Differential measurements by gating on the A — 1 
residual nuclei which will provide information on the annihilation of the incoming p on a single target neutron or 
proton, respectively, are being explored [TBI ]. Parallel to the nuclear structure studies the experimental program 
should also envision investigations of elastic pA scattering. A consistent and much extended set of clastic angular 
distributions is highly desired by theory for a better understanding of antiproton interactions in nuclear matter. 

Supported in part by DFG, contract Le 439/05, GSI and BMBF. 
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